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Abstract 



First principles FLMTO-GGA electronic structure calculations of the new 
medium- Tc superconductor (MTSC) MgB2 and related diborides indicate that su- 
perconductivity in these compounds is related to the the existence of p^^y-band 
holes at the F point. Based on these calculations, we explain the absence of 
medium- Tc superconductivity for BeB2, AIB2 ScB2 and 1^-62- Fhe simulation of a 
number of MgB2-hased ternary systems using a supercell approach demonstrates 
that (i) the electron doping of MgB2 (i.e., MgB2-yXy with X = Be,C, N,0) 
and the creation of isoelectronic defects in the boron sublattice (nonstoichiometric 
MgBy^2) are not favorable for superconductivity, and (ii) a possible way of search- 
ing for similar MTSC should be via hole doping of MgB2 (i.e., Mgi-.xMxB2 with 
M = Be, Ca, Li, Na, Cu, Zn) or CaB2 or via creating layered superstructures of 
the MgB2/CaB2 type. A recent report of superconductivity in Cu doped MgB2 
supports this view. 
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The discovery of superconductivity in magnesium diboride {MgB2)^ has attracted a great 
deal of interest in this system. The transition temperature for MgB2 {Tc ~ 39 K) exceeds 
by almost two times the record values of Tc for conventional Bl- and AlS-type intermetallic 
superconductors (SC)i. As distinct from the high-temperature SC, MgB2 has an exclusively 
simple composition and crystal structure!. 

The finding of a new medium-Tc SC (MTSC) raises a number of questions: 1. What is 
the nature of the superconductivity in MgB2'? 2. Is MgB2 a unique compound of this kind 
or is it the first representative of a new class of MTSCs? 3. What physical and chemical 
properties should related MTSCs possess? 

It is worth noting that among boron-containing phases MgB2 is not a unique supercon- 
ductor. A systematic search for superconductivity in a wide range of d-metal borides (and 
other metal-like compounds of a transition metals (TM) with light non-metals (carbides and 
nitride^l )) has been carried out earlieriH. For instance, it was shown that the supercon- 
ducting transition temperature for TM diborides {MB2, M = Ti, Zr, Hf, V, Ta, Cr, Mo) is 
below ~ 0.4 K. Only NbB2 was found to be a superconductor with a Tc of about 0.6 K. 
A large class of superconducting (with Tc of several K) ternary {LnRuB2, LnRhiB^) and 
pseudo-ternary (Lni^^Lny) Rh^B^ borides is known0. In 1994, Cava et ali discovered super- 
conductivity {Tc ~ 16 — 23 K) in a new class of intermetallic borocarbides (IBC), namely 
layered compounds of the LnM2B2C type. Numerous theoretical and experimental studies 
(see reviewi ) made it possible to place IBC into the class of conventional BCS superconduc- 
tors. It is significant that bands near the Fermi level of IBC participating in electron-phonon 
interactions (with a high-frequency Baig mode) are determined by interatomic a-bonds in 
the TMB4 tetrahedra. For all the above mentioned borides, the most important role in 
their superconductivity is played by d-electron atoms. 

The electronic properties of MgB2 are of quite a different kind. Earlier band structure 
calculations of MgB2 (by semiempirical LCAO0 and first principles full potential linear 
muffin-tin orbitalS, FLMTO, methods) showed that the upper filled energy bands of MgB2 
are formed mainly due to strong B — B interactions (in honeycomb layers of boron atoms). 



Analogous conclusions were made recently by other authorslllhEl, who pointed out that the 
coexistence of 2D in-plane and 3D interlayer bands is a peculiar feature of MgB2- Based 
on an estimate of phonon frequences and band structure calculations, Kortus et al0 explain 
the superconductivity in MgB2 as a result of strong electron-phonon coupling and An and 
Pickett0 attributed it to the behavior of p^^^^-band holes in negatively charged boron planes. 
All authors0^ii emphasize the most significant role of metallic B states in the appearance of 
superconductivity. According to the McMillan formula for Tc-B, the most probable MgB2- 
based MTSC should have a high density of states at the Fermi level, N{Ef), high averaged 
electron-ion matrix elements, as well as high phonon frequences, which increase for light 
elements and depend on B — B and M — B bonding. 

To examine the possibility of superconductivity in related diborides, we studied the band 
structure of MgB2 in comparison with that of CaB2, BeB2, AIB2, ScB2 and YB2 which 
are typical representatives of different groups of AlB2-like diborides formed by s-, p-, and 
(i-metals, respectively. This makes it possible to analyze the following effects in: (1) MgB2 
versus CaB2 and BeB2, which are isoelectronic and isostructural. The basic changes will 
be due to structural factors (changes in lattice parameter, a, and interatomic distance, c/a). 
(2) MgB2 vs. AIB2, for which the main differences are expected as a result of changes in the 
filling of bands with addition of an electron. (3) MgB2 vs. ScB2 and YB2, where alongside 
changes in the filling of bands, the band structure and interatomic bonding will depend on 
the M-sublattice (s- or d-metal). (4) AIB2 and ScB2, YB2, which are also isoelectronic and 
isostructural, and changes will be due to the nature of the M-sublattice {p- or d-metal). 

In addition, using the supercell approach, we have carried out a theoretical search for 
possible superconductors among some M(7i?2-based ternary systems. For this purpose, we 
modelled the effect produced on the band structure of MgB2 by (1) boron sublattice doping 
(with Be, C, N, O impurities), (2) magnesium sublattice doping (with Be, Ca, Li, Na, Cu, 
Zn impurities), and (3) the presence of lattice vacancies in Mg- and -B-sublattices. 

The above diborides have a hexagonal crystal structure (yl/_B2-type, space group 
PQ/mmm, Z = 1) with atomic positions M (a) 000, B (d) 1/3,2/3,1/21. The lattice 



is composed of layers of trigonal prisms of M-atoms in the center of boron atoms, which 
form planar graphite-like networks. Their band structures were calculated by the FLMTO 
method within the generalized gradient approximation for the exchange -correlation po- 
tential (FLMTO-GGA)0. Ternary and nonstoichiometric diborides were simulated with 
12-atom supercells (2x2x1) and experimental lattice parameters for MgB2, BeB2, AIB2, 
ScB2 and YB2 taken from Ref.i. For the hypothetical phase of CaB2, the lattice constants 
(a = 3.205 A, c/a = 1.24) were determined by total energy minimization. 

MgB2. The band structure and DOS are shown in Figs. |l] and |^. The high energy part 
of the valence band (VB) of MgB2, made up predominantly of B 2p-states, form two distinct 
sets of bands of the a {2px^y)- and vr {pz)-tjpes, whose k dependence differs considerably. 
For B 2px,y the most pronounced dispersion is observed along V — K. These bands are of the 
quasi-two dimensional (2D) type, form a fiat zone in the direction (F — A) and reflect the 
distribution of ppo--states in the B layers. These states make a considerable contribution to 
N{Ef), forming metallic properties of the diboride with Ep located in the region of bonding 
states, the conductivity is due to hole carriers. 

The B 2p2-like bands are responsible for weaker pp^-interactions. These 3D-like bands 
have maximum dispersion along {T — A); Mg s,p- and B s-states are admixed with B 2p 
bands near the bottom of the VB and in the conduction band. Thus, the peculiarities of 
the electronic properties of MgB2 are associated with metallic 2p-states of B atoms located 
in planar nets, which determine the DOS in the vicinity of Ep. 

The valence charge density (VCD) maps (Fig. ^) demonstrate that the interatomic bonds 
in the boron graphite-like layers are highly covalent, the Mg-Mg bonds have metallic char- 
acter and then there are the weak interlayer covalent Mg — B interactions. Estimates of in- 
teratomic bonds strengths (with FLMTO calculations using empty Mg and B sublattices0) 
also shows the largest contribution to the cohesion energy (Ecoh) from B — B interactions 
(B-B (68%), B - Mg (23%) and Mg - Mg (9%)). In view of these resuhs (see also&El), 
let us discuss the possibility of superconducting properties in related diborides. 

MgB2 versus CaB2. The electronic bands of these isoelectronic (electronic concentration 
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per atom, Ue = 2.67) and isostructural compounds turned out to be similar on the whole, 
see Fig. 0. The differences are due to the increased lattice parameter, a, and interlayer 
[M — B) distances determined by M ionic radii [vc = 0.74{Mg) and 1.04 A (Ca)). The main 
change is a downward shift of the p^-bands as compared to px,y-ha.nds, so that the crossing 
point is located at Ep (Fig. |ip and N{Ef) increases to the highest value of any material 
given in Table |l|. This may well indicate that a possibly higher Tc may be obtained if CaB2 
were stabilized. To this end, we also studied the lattice stability of the hypothetical CaB2 
by calculating the formation energy as a difference in the total energy with reference to 
the constituent elements in their stable modifications, viz. hep calcium and rhombohedral 
boron {a — B12). It was found that the formation energy has a small, but negative value 
{Eform = —0.12 eV/f.u.), which may be indicative of the possibility of real synthesis of 
Cai?2- For comparison, Eform for MgB2 was found to be -1.21 eV/f.u.. 

Let us outline ways to stabilize CaB2. It is known for stable diborides that c/a does 
not exceed ~ 1.165^. For the equilibrium state of CaB2, we obtained c/a = 1.24. The most 
obvious way of decreasing the interlayer distance consists in the partial replacement of Ca 
by atoms with smaller radii. To obtain superconductivity in such doped ternary systems, 
these metals should be either isoelectronic with Ca (for example, Mg, Be) or hole dopants 
{Li, Na, etc.) — as discussed below. They may also be prepared as layered superstructures, 
for example, .../ Ca/ B2/ Mg / B2/ .... 

MgB2 versus BeB2. Their band structures turned out to be similar, see Fig. The 
differences are due to a downward shift of p^^^^-bands, so that they are absent above Ep 
at r. This brings about changes in the Fermi surface topology: cylinders along T — A are 
transformed into cones. As a result, BeB2 is not a MTSC — as supported by the absence 
of superconductivity in recent experiments^. 

MgB2 versus AIB2 . For AIB2 {ug = 3.00), the ^^..^^-bands are completely filled and 
AIB2 is not a superconductor. A recent first report^ on the electron-doped Mgi-yAlyB2, 
demonstrates that Tc decreases smoothly and vanishes a.t y = 0.25. According to our 
FLMTO estimates^, the interatomic bonding exhibits some changes: the contribution of 



Al — Al bonds to the cohesive energy increases by almost two times and the M — B bonds 
become covalent (see AIB2 VCD maps in Fig. 4 of Ref.§). 

MgB2 versus ScB2 and YB2. Quite a number of studies of the electronic properties of 
these diborides are known to date@^^. We previously0@ performed FLMTO calculations of 
all hexagonal diborides of 3d- {Sc, Ti...Fe), Ad— {Y, Zr...Ru), and 5d— {La, H f...Os) metals 
and analyzed variations in their chemical stability and some other properties (e.g., melting 
temperatures, enthalpies of formation). We found that the evolution of their band structures 
can be described within a rigid-band model (RBM). For M = Ti, Zr, Hf {rie = 3.0), Ep 
falls near the DOS minimum (pseudogap) between the fully occupied bonding bands and 
unoccupied antibonding bands. In the RBM, MgB2 and Sc, Y, La diborides all have 
partially unoccupied bonding bands. Comparison of their band structures shows that for 
ScB2, the 2D-2px y bands are almost filled and the hole concentration is very small (near 
A in Fig. H). The Sc rf-band along F — M is below Ep and the large contribution to DOS 
at Ep {N{Ep)=1.06 states/eV-f.u.) is due to Sc ci-states. The covalent M — B bonding 
increases considerably and the B 2pz-\ike bands are shifted downwards at K. For both ScB2 
and YB2 there is a small hole concentration of B 2px,y states at A. Thus, one can expect 
for these diborides only low temperature superconductivity (for ScB2 Tq ~ 1.5 kII). 

Modeling of new MgB2-based ternary borides. In searching for possible superconductors 
among the M(7i?2-based ternary systems, we simulated the effect produced on the band 
structure of MgB2 by (i) doping the B sub lattice with C, N, O impurities, (ii) doping 
the Mg sublattice with Be, Ca, Li, Na, Cu, Zn impurities, and (Hi) adding vacancies in 
the Mg— and B— sublattices — all of which can be divided into three groups: electron-, 
hole-dopants and isoelectronic defects. 

The above mentioned RBM is widely used to analyze such substitutions. Based on the 
DOS of MgB2 (Fig. 0), this model determines the following effects: hole doping should lead 
to partial "depopulation" of bonding bands, a shift of Ep deep into the VB, and an increase 
in N{Ep). A reverse effect, namely a Fermi level shift to the region of the DOS minimum 
(pseudogap between bonding and antibonding states), can be expected for electron doping. 



The role of isoelectronic dopants remains unclear within this model. One can only suppose 
that in this case the main effects will be associated with lattice deformations, i.e., with the 
changes in a and c/a. By contrast, our FLMTO calculations for hypothetical ternary and 
nonstoichiometric borides demonstrate a more complicated picture (see Figs. ^ and ^ and 
Table |I|) and lead us to the following conclusions: 

Boron suhlattice doping. As Ue decreases {MgB2 MgBu^), Ep is shifted to higher 
binding energies, but the value of N{Ef) becomes smaller due to a partial breakdown of -B — 
B bonds in the graphite-like B layers and changes in the energy band dispersion near Ep- For 
MgBij'^Co_25, all bonding states are completely filled and Ep falls in the pseudogap, hence a 
decrease in N{Ef). A further increase in {MgBi Y^{N, 0)0.25) results in the occupation of 
2p-antibonding states. Thus, we do not see any prospects for obtaining superconductivity in 
M(yfi?2-based ternary compounds by doping the i?-sublattice. The considerable energy cost 
of breaking the strong B — B bonds should be emphasized. As reported in0, the contribution 
oi B — B bonds to the cohesive energy of MgB2 is 3.43 eV/atom. The corresponding values 
for Mg — Mg and Mg — B bonds are 0.45 and 1.16 eV/atom, respectively. This means that 
substitutions in the Mg sublattice will be energetically more preferable. 

Magnesium sublattice doping. As decreases in the series MgB2 — >■ 

Mgo,'r5Li{Na)o,25B2 MgQ,-j^B2, MgQ,^Li{Na)o,^B2, Ep shifts to higher binding energies, 
but no increase in N{Ep) is observed — as predicted within the RBM, see Table |I]. Thus, 
we found that for Li and Na substitution, N{Ep) is almost independent of the hole dopant 
concentration. This fact is in agreement with the results of An and PickettS. The depen- 
dence of the energy of (T(2pj. y)-bands (forming cylinders along T — A) at T relative to Ep 
{E'^{r — Ep)) is shown in Fig. ^. This dependence for both the hole and electron doping of 
MgB2 is close to linear. It is seen that hole doping in M(yfi?2-based ternary systems serves the 
peculiarities of band structure of MgB2 near the Ep. Moreover, both N{Ep) and E'^{r~Ep) 
in these systems are practically independent of the hole dopant type, and are determined 
only by the hole concentration (see the results for MgQ -j^LiQ 25B2 and MgQ j^NaQ 25B2 or 
for M^fQ. 5^20.5-82, MgQ,^NaQ,^B2 and M(7o.75i?2 in Table |l] and Fig. ||). For the isoelectronic 
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substitution MgB2 Mgo_j5Be{Zn)o,25B2, the value of N{Ef) decreases but pa;,y-bands 
move up by 0.1 eV at F comparing with MgB2, Fig. ^. 

Very interestingly, in a recent experiment^ Tq was found to be increased up to 49 K with 
nominal 20 % Cu doping of MgB2. In our calculation, we found that 25 % Cu substitution 
for Mg leads to an increase of N{Ep) (c.f.. Table |I]) and here the top of p^^^^-bands at F has 
the same location as for other doped systems with the same hole concentration (c.f.. Fig. 

This adds strong support to the validity of the systemetic results presented here. 

We believe that further searches of new MTSC should be made via hole doping of MgB2 
(or BeB2, CaB2) or creating layer superstructures of the MgB2/CaB2 type. We are contin- 
uing these investigations taking into account the crystal lattice relaxation and calculating 
the substitution energies for impurities in ternary systems. 

We thank V. Antropov, J. Jorgensen and D. Novikov for useful discussions. Work at 
Northwestern University supported by the U.S. Department of Energy (Grant No. DE- 
F602-88ER45372). 
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FIGURES 

Fig. 1. Band structures of (a) MgB2, (b) BeB2, (c) ScB2, (d) CaB2, (e) AIB2, and (f) YB2. 



Fig. 2. Density of states MgB2, CaB2, Mg'0.75 -^00.25-^2 and -^50.75^^^0.25-62- 

Fig. 3. The valence charge density (VCD) maps for MgB2 in in-plane (B — B, Mg — Mg) and 
interlayer {Mg — B) sections. 

Fig. 4. The energy of the a{2px,y)-han.d at F relative to Ep {E'^{T — Ep)) as a function of 
the hole- or electron-type doping {n^ - the difference of valence electron concentration relative to 
MgB2). 
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TABLES 

Table 1. Density of states at the Fermi level {N{Ef), in states/(eV-f.u.)) for MgB2 and some 
related borides. 



System 


N(Ef) 


System 


N{Ef) 


MgB2 


0.73 


Mc/o.75^«0.25-B2 


0.73 


CaB2 


0.92 


Mgo.5Lio.5B2 


0.75 


BeB2 


0.47 


MgQ.75NaQ.25B2 


0.73 


MgBi,Y5Co.25 


0.56 


MgQ.5NaQ.5B2 


0.76 


MgBi,-r5 


0.61 


Mgo.75cuo.25B2 


0.89 


Mgo,75B2 


0.74 


MgQ.75BeQ.25B2 


0.62 






MgQ.75znQ.25B2 


0.51 
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